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Abstract 

Aims. Generally, the oblateness of a planet or moon is what causes rings to settle into its equatorial plane. However, the 
recent suggestion that a ring system might exist (or have existed) about Rhea, a moon whose shape includes a strong 
prolate component pointed toward Saturn, raises the question of whether rings around a triaxial primary can be stable. 
We study the role of prolateness in the behavior of rings around Rhea and extend our results to similar problems such 
as possible rings around exoplanets. 

Methods. Using a Hamiltonian approach, we point out that the dynamical behavior of ring particles is governed by three 
different time scales : the orbital period of the particles, the rotation period of the primary, and the precession period 
of the particles’ orbital plane. In the case of Rhea, two of these are well separated from the third, allowing us to average 
the Hamiltonian twice. To study the case of slow rotation of the primary, we also carry out numerical simulations of a 
thin disk of particles undergoing secular effects and damping. 

Results. In the case of Rhea, the averaging reduces the Hamiltonian to an oblate potential, under which rings would 
be stable only in the equatorial plane. This is not the case for lapetus; rather, it is the lack of a prolate component to 
its shape that allows lapetus to host rings. Plausible exoplanets should mostly be in the same regime as Rhea, though 
other outcomes are possible. The numerical simulations indicate that, even when the double averaging is irrelevant, 
rings settle in the equatorial plane on an approximately constant time scale. 
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1. Introduction 


The oblateness of a planet (i.e., the slight flattening 
along its spin axis), characterized by the coefficient J 2 , de¬ 
fines t he Laplace plane in w hich potential rings may be 


stable (Tremaine et al. 2009). However, some bodies, like 


moons orbiting their planet, are also elongated by tidal ef¬ 
fects in the direction toward the object they are orbiting. 
This elongation is described by the coefficient (722- Its ef¬ 
fect on the stability of planetary rings is the subject of this 
work. 


In 2008, the Cassini spacecraft observed some sharp 
symm etric drops in the flow of charged particles around 


Rhea (Jones et al.l 2008). This observation was interpreted 


as being due to the presence of narrow rings orb iting the 
moon. Eurthermore, Schenk and McKinnon (2009) found a 
narrow equatorial color signature that was ascribed to the 
fallout of ring material and was cited a s furth er evidence of 
Rhean rings. However, Tiscareno et al. (2010) ruled out the 
current existence of rings with an in- depth imaging study 
of the equatorial region, after which Schenk et al. (2011) 
suggested that the narrow equatorial color signature may 
still indicate the existence of rings in the past. 


Rhea would have been the first triaxial object known 
to possess a ring system. As for most moons close to their 
planet, both Rhea’s J 2 and C 22 are significant. This work 
is motivated by the question of whether it is even possible 
for rings to be stable in such an environment. 


The calculations we develop can be adapted to lapetus, 
for which rings have been invoke d to explain the origin of its 


equatorial mountain ra nge (Ip 2006 


Levison et al. 


2011 


Bombard et al. 2012). lapetus is different from Rhea in 
that it is simply oblate, but a comparison with our model 
remains interesting. 

Some exoplanets, such as hot Jupiters, are likely to be 
drawn out by tidal effects in the direction of their star. This 
deformation, added to the fast rotation of the planet, must 
have complex effects on the presence of putative rings. 

We will carry ou t a Hamiltonian analysis, following the 
model devel oped by [Dobrovolskis et al. (|1989) for N eptune 
polar rings. Steiman-Cameron and Durisen (1984) used a 
similar model to describe gas rings surrounding triaxial ga¬ 
laxies, which is better adapted to the current case of slowly 
rotating bodies. 

In the next section, we show that the inertial frame is 
better adapted to the description of the motion than the 
rotating frame fixed on Rhea. In Sect. 3, we derive the 
perturbative Hamiltonian and show that the prolateness 
has no average effect. In Sect. 4, we investigate the case of 
lapetus and exoplanets. Section 5 contains a computational 
simulation that allows us to study slower rotations. 


2. The choice of a frame 

We consider a test particle moving around Rhea on a 
nearly-Keplerian orbit perturbed by the shape of Rhea and 
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the presence of Saturn. In this work, the particle’s orbit is 
assumed to be circular. This is probably a good approxi¬ 
mation for any putative rings. 


2ti I dTT^a^ 


2.1. Definition of the frame and time scales 

We will first define a system of axes in which the expres¬ 
sions for the perturbations are simple (Fig.[^. In general, I 
is inclination, D is longitude of ascending node, and 0 and 
(j) are the particle’s colatitude and longitude, respectively. 
The xr axis always points toward Saturn during the mo¬ 
tion of Rhea; therefore, it corresponds to Rhea’s long axis 
and to the C 22 term. The zr axis is taken parallel to Rhea’s 
angular momentum, i.e., perpendicular to its orbital plane; 
it is the short axis. The ^r axis completes the other two to 
form a direct basis. The angles in this frame will be labeled 
with a subscript R : /r, f^R, etc. 


ZR 



( 1 ) 


GMr \ r ^) 


where a is the semi-major axis and n = is the 

mean motion for the test particle on its orbit around Rhea, 
Ur and riR are the same for Rhea on its orbit around Saturn, 
Mr and Ms are the respective masses of Rhea and Saturn, 
Rr is the radius and J 2 is the oblateness parameter for 
Rhea, and G is Newton’s constant. 

In Tp, we esti mated Ct as O = 3/2n J 2 ( Rr/u)^ from 


Equation 6.250 in Murray and Dermott ( 1999| ). This is valid 


when the only perturbation is oblateness, but it should give 
the correct precession rate to the first order in our situation, 
as the J 2 and C 22 are almost equally significant. 

In the case of Rhea, J 2 = 889 10-^ Mr = 2.308 10^^ 





Schubertl 

2007). If we take a 1.5 Rr, as in | Jones et al. 

wk. 


J yR 

to 

0 

0 

00 

8 time scales are as given in Table 


Table 1. Time scales of the problem for Rhea at a = 1.5 Rr 


Figure 1. Angles in the rotating frame 


T 

0.227 d 

Tr 

4.52 d 


384 d 


This frame is well adapted to express Rhea’s gravitatio¬ 
nal potential, but it will not allow a simple description of 
the behavior of putative rings. In order to better unders¬ 
tand the importance of the frame, we have to focus on the 
three relevant time scales of our problem : r, the period of 
one orbit of the test-particle around the moon; tr, the per¬ 
iod of one orbit of the moon around the planet (identical to 
the rotation period of the moon); Tp, the precession period 
of the test particle’s orbit due to the non-sphericity of the 
moon. This precession is responsible for the settlement of 
the rings in preferred planes, so we can take this time as a 
low estimation of the settling time. 

If Tr ^ r and tr ^ Tp, the particle has time to settle in 
the preferred planes relative to the moon, and the rotating 
frame coordinates is the most convenient. The rings turn 
with the moon during its motion. 

If Tp ^ Tr, the rings remain steady in the inertial frame 
of Saturn during one period of Rhea. Therefore, on the time 
scale on which the particles can settle into rings, they only 
see an average potential : the rotating C 22 term averages 
to give an azimuthally symmetric shape to the satellite. In 
this case, the frame that rotates with Rhea is not the most 
convenient. 


2.2. The case of Rhea 

We estimate the different time scales as 


Clearly, Tp ^ tr. Thus, the most convenient frame is 
the inertial frame centered on Rhea that does not rotate. 
The transformation between the inertial and the rotational 
frames is just a rotation of angle n^t around the 2 ;-axis. 

There is a factor of at least 50 between each time scale. 
This invites us to make a double averaging : the first on the 
orbit of the particle around Rhea, and the second on the 
orbit of Rhea around Saturn. 


3. The perturbative Hamiltonian 

3.1. Derivation of the Hamiltonian 


The system is conservative, so the semi-major axis a is 
constant. Therefore, an orbit is described by its inclination 
/r and the position of the ascending node ffp. Another set 
of coordinates is relevant here, the particle’s colatitude and 
longitude (^r,0r), because the perturbative Hamiltonian 
is expressed in thei r terms. The time-evolution of /r and 
f^R is governed by (Kaula 


1968) 


= 




1 


dH 


na^ sin(/R) 

1 dH 


( 2 ) 


na2 sin(/R) din 

All Hamiltonians will be written per unit of the test 
particle’s mass. The most important perturbation for us 
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is the one that results from the shape of Rhea (Steiman- 


Cameron and Durisen 1984) 


= n^Rl 


J 2 


2 


“ o “ 3C22 sin^ 6 >r cos 20 r 


( 3 ) 

The particle also feels the effect of the gravitational field 
from Saturn. The quadrupole term (i.e., order 2 in the de¬ 
velopment in a/ap) is a good approximation for this per¬ 


turbation. It is given by (Murray and Dermott 1999) 


-ffs = -3 sin^6»RCOs^(/)R)-nlal. (4) 

We will now take the average of the perturbation over 
the duration r = 27r/n of one orbit of the particle around 
Rhea, H = (1/r) fj Hdt, where H = i^R + Hs. Dropping 
the parts of the Hamiltonian that depend neither on I nor 
on D, we obtain 


^ n^Rl [3J2 . 


H 


sim I — 3 C 22 cos 2(D — npt) sim I 


+ 


3n‘^a‘^ 


sin^(D — npt) sin^ /. (5) 


We employ a similar method to average the Hamiltonian 
once more, over a time tr, which corresponds to one orbit of 
Rhea. Now the C 22 term vanishes, and the double-averaged 
Hamiltonian is 






sin^ I. 


(6) 


where q is the rotation parameter defined by 


V Rr / 


( 9 ) 


Figure [^represents the space (J 2 , ^/q). The gray region 
corresponds to the region where Eqs. are satisfied, and 
the dark gray zone to the sets of parameters for which each 




Figure 2. Phase space in terms of Eqs. Rhea lies in the 
dark gray zone, which means that we can safely apply our 
double-averaging method. lapetus is not in that region. The 
dashed curve gives the position of homogenous, synchro¬ 
nous ellipsoidal objects. 


3.2. Interpretation for Rhea 

The sin^ / factor is characteristic of a J 2 perturbation. 
The presence of Saturn is equivalent to an increase in the 
oblateness of Rhea, and can be taken into account with an 
effective J 2 . Erom Eq. this effective J 2 is given by 


Jt” = J2 + 




( 7 ) 


Close to the surface of Rhea, i.e., for a = Rr, the second 
part is already equal to 375 10“^. This means that the effect 
of Saturn quickly becomes dominating when a increases. 


4. Extension of the model 


4.1. Validity of our approximations 


The averaging over several orbits of the particle around 
Rhea requires that tr be much larger than r. The other 
averaging, over several revolutions of Rhea necessitates that 
Tp be much larger than tr. We will study the most favorable 
possible case of a particle close to the moon’s surface : a ^ 
Rr. With the definitions given in Eq. [^we can rewrite the 
two previous conditions as 


^/q I 




< 1 ’ 


( 8 ) 


4.2. lapetus 


lapetus’ J 2 has never been measured, but Thomas et al. 


(20071) give a = 747.4 km, h = c = 712.4 km, and mean ra ¬ 
dius Ri = 735.6 km, yielding J 2 = 1890 10“^ (Yoder 1995), 
under the simplifying assumption that lapetus is homoge- 
nous. With the definitions of the time scales given in Eq. 
we calculate the values of Table at a = 1.5 Ri. 


Table 2. Different time scales for lapetus at a = 1.5Ri 


T 

0.248 d 

n 

79.3 d 


197 d 


The condition required to make the averaging over seve¬ 
ral lapetus rotations is not well satisfied : n/r^ = 0.4 (see 
Eig.l^. This would be an important effect if lapetus did not 
have an azimuthal asymmetry. Thus, equatorial rings can 
be stable around Rhea and lapetus, but for different rea¬ 
sons : for lapetus it is the absence of a significant C 22 term, 
while for Rhea the C 22 term is neutralized by the fast re¬ 
volution. If lapetus’ C 22 had the same magnitude as its Jg , 


polar rings would also be stable (Dobrovolskis et al. 1989|), 
because the moon’s rotation is slower than the precession 
of a particle in orbit around it. 
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4.3. Exoplanets 


We now consider whether the non-spherical shape of 
some exoplanets (such as hot Jupiters) can destabilize a pu¬ 
tative ring system. In such an environment, any ring system 
would have to be composed of refractory material, not ice 
(Schlichting and Chang 2011). We choose a simple model 
of a homogenous synchrono us planet at hydrostatic equ ili- 


brium. Therefore 


we have (Murray and Dermott 1999| ) 


J2 = (10) 

This relation further simplifies the problem : the condi¬ 
tion for averaging the exoplanet’s shape over its rotation 
(^p/^R ^ 1) then becomes ^/q <C 3, which is automatically 
satisfied with the first condition for averaging the particle 
along its orbit (tr/t 1 or equivalently y/g 1). We re¬ 
present the hydrostatic relation graphically in Fig. with 
the dashed curve = qj2. 

We express q in terms of more physical parameters : the 
mean density of the exoplanet the mass of the central 
star M*, and the semi-major axis of the exoplanet a^. If we 
want a factor of 10 at least between each time scale, we can 
look for the limit parameter when the two others are fixed. 
The restrictive condition can then be written as 


^ TT 

75’ 


( 11 ) 


If is set, there is a minimal mean density for the 

exoplanet below which we cannot apply our averagings : 
^min ^ 75M*/7rag. Figure I shows the minimum density 
of the exoplanet as a function of the semi-major axis for 
several given masses of the star from 0.1 Mq to 10 Mq. 


Pe™ (10^ kg/m^) 



Figure 3. The minimum density for which double¬ 
averaging is appropriate, as a function of the semi-major 
axis. Contours are for different masses of the central star. 


Figure indicates that above approximately 0.1 AU, 
all exoplanets are dense enough to satisfy the conditions. 
For = Mq and a between 0.012 AU and 0.022 AU, 
only the planets with a density greater than Jupiter’s are 
dense enough. For even smaller values of the semi-major 
axis, only telluric planets have the required density. Very 
few exoplanets orbit so close to their star that no reasonable 
density can fit the criteria. For hot Jupiters, the typical 


values are about Ue = 0.05 AU and = Mq, so the 
minimum density equals 113 kg/m^. This value is 

even lower than that of HD 209458 b, which is one of the 
l east dense of the kn own exoplanets, with ^ = 350 kg/m^ 
(Brown et al. 2001). For these exoplanets the conclusions 
about the role of the prolateness are the same as for Rhea. 


4.4. Hill's sphere limit 

In the case of exoplanets, the star also directly attracts 
the test particle. Hill’s sphere radius can be estimated by 
Rh = ae(Me/3M*)^/^. It would have no meaning to consi¬ 
der rings outside Hill’s sphere. Therefore, in a system that 
is able to shelter rings. Hill’s radius must be greater than 
a, the semi-major axis of the orbiting particle. Setting 
a = aRe, this condition can also be written in terms of 
Q : 


q< 


1 


( 12 ) 


Even for very close rings, the maximum allowed value 
of q is on the order of 0.1. This is the most unfavorable 
case, and yet it nearly already satisfies our hypotheses for 
the double averaging (y/g ^ 1). To sum up, we should be 
careful when g increases too much, but in these cases the 
attraction of the central body is likely to prevent any rings 
from forming. 


5. Computational simulation for slower rotations 

5.1. Model and results 

Not all planets are homogenous or at hydrostatic equili¬ 
brium. lapetus is a typical example of a planet whose J 2 is 
too large to be caused by rotation at its distance to Saturn. 
Thomas et al. (2007) called this a “fossil bulge”, based on 
the likely idea that it is left over from a time when lapetus 
rotated faster. Therefore we consider what happens when 
the rotational period of a moon (e.g., Rhea) tr is not short 
compared with the precession time Tp. The Hamiltonian to 
consider is then the one given by Eq. Given its com¬ 
plexity, we will use a simulation. 

Our simulation will describe the behavior of a thin disk 
of interacting particles. To take into account these inter¬ 
actions, we will s imply give onr part i cle di sk a kinematic 
viscosity u, as in Dobrovolskis et ^ ( 1989| ) and Steiman- 
Cameron and Durisen ( |1988p . We modify the system of Eqs. 
2 by adding a dissipative te rm (/d,^d) (see Eqs. 9 and 10 
in Dobrovolskis et al. 1989) : 


= 


Un = — 


a'^I ^ r/'“V 


mdl 

^ +2cotJ—— 


(13) 


We assume that a is constant. Indeed, the settling times 
are much shorter than the amount of time where the inflow 


is significant (Dobrovolskis et al. 1989). The ring is then 
modeled as a set of 100 annular elements. Eorty kilometers 
separate the innermost annulus from the outermost, and 
all annuli have the same initial inclination. Then, we take 
a time step of tr/IOO. The kinematic viscosity is estimated 
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to be about 1 m^/s , using the model of Steiman-Cameron 
and Durisen (1988), with ice particles of millimeter size. 
The settling time is denoted Tg. 

Testing our simulation with an object having the cha¬ 
racteristics of Rhea, we find satisfactory results. No matter 
what the initial inclination is, the rings always settle into 
the equatorial plane, as expected. The settling time we ad¬ 
mitted for our theoretical study was equal to the precession 
period Tp. The simulation gives settling times Tg which are 
closer to ten times Tp. However, our goal was only to ob¬ 
tain a theoretical behavior for the rings, not to predict the 
detailed settling process. Such a difference is therefore ac¬ 
ceptable. 

We can now try to vary the parameters of the problem, 
and especially to study what happens when tr Tp. We 
will decrease the rotational frequency np. This is equivalent 
to moving Rhea outward from Saturn while retaining its 
current shape (J 2 and C 22 ) as a fossil bulge. We will call this 
new frequency to avoid confusion with Rhea’s proper 
characteristics. 

We find that the behavior of the rings is roughly inde¬ 
pendent of n™. The rings still settle in the equatorial plane 
during a nearly constant duration Tg. Table shows that Tg 
is not dependent on the rotational speed. 


to shelter rings owing to their small Hill spheres. Therefore, 
results obtained for Rhea can be extended to almost all syn¬ 
chronous exoplanets. 

Finally, we carried out a computational simulation focu¬ 
sed on what happens when the rotational period becomes 
as significant as the precession time. This happens when 
the J 2 is larger than it should be, as it is in the case of 
lapetus. Though it neglects all effects that are not azimu- 
thally averaged, our simulation indicates that, over a wide 
range of rotational periods, a thin particle disk will become 
rings in the equatorial plane. The duration for the settling 
is roughly independent of the rotational period. 

Thus, in all cases that were studied, the prolate shape 
seems to have no determining effect on the shape of the 
rings, keeping in mind that classical shape models for the 
exoplanets were assumed. Nothing prevents the discovery 
of more original shapes, in particular higher prolate 
deformations allowing more imaginative ring systems. 

We thank Joe Burns, Phil Nicholson, Matt Hedman, 
and Tom Steiman-Cameron for helpful conversations. We 
acknowledge funding from the Ecole Normale Superieure 
de Cachan and from NASA Outer Planets Research 
(NNX10AP94G). 


Table 3. Settling time at several rotational frequencies. 


Ts (10^ days) 

nT/nn 

5.6 

10“^ 

5.3 

10“^ 

5.3 

1 

5.4 

10 

5.0 

10^ 


5.2. Caveat 

The simulation results should be interpreted with some 
caution. While the simulation robustly reproduces secular 
effects, it may not represent any effects that are not azimu- 
thally averaged, such as any imparting of eccentricity to the 
ring or any resonant effects due to the rotating character 
of the C 22 term. 


6. Discussion 


The separation between the different time scales in the 
case of Rhea allows us to carry out a double averaging and 
obtain a more understandable Hamiltonian. We have shown 
that rings can settle only in the equatorial plane. This rules 
out the possibility of rings in a region around Rhea that 
would not have been carefully investigated. The results of 
Tiscareno et al. (2010|) then seem to definitively exclude the 
present existence of tnese rings. 

Our model turned out to be interesting for other celes¬ 
tial objects. The most likely objects to have a shape close to 
Rhea’s are exoplanets orbiting close to their stars. Nearly 
all of them should be suitable for our model, except for 
those orbiting at less than 0.02 AU from their star. A weak 
density could put such rare planets into a more complex re¬ 
gime. On the other hand, such close exoplanets are unlikely 
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